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Abstract 

We introduce the notion of relative volatility/intermittency and demonstrate how 
relative volatility statistics can be used to estimate consistently the temporal vari- 






ation of volatility/intermittency even when the data of interest are generated by a 
non-semimartingale, or a Brownian semistationary process in particular. While this 
estimation method is motivated by the assessment of relative energy dissipation in 
empirical data of turbulence, we apply it also to energy price data. Moreover, we 
develop a probabilistic asymptotic theory for relative power variations of Brownian 

i_i semistationary processes and ltd semimartingales and discuss how it can be used for 

f-H inference on relative volatility/intermittency. 
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1 Introduction 

The concept of (stochastic) Volatility is of central importance in many fields of science. 
In some of these the term Intermittency is used instead of volatility. Thus volatil- 
ity/intermittency has a central role in mathematical finance and financial econometrics 
[BNS10], in turbulence, rain and cloud studies [LS06, Way06] and other aspects of environ- 
mental science [PA12], in relation to nanoscale emitters [FVKJ13], magnetohydrodynamics 
[MP09] , and to liquid mixtures of chemicals [Sre04] , and last but not least in the physics 
of fusion plasmas [CSVMAVOO]. In turbulence the key concept of Energy Dissipation is 
subsumed under that of intermittency 1 . 

Speaking generally, volatility/intermittency is taken to mean that the phenomenon 
under study exhibits more variation than expected; that is, more than the most basic type 
of random influence 2 envisaged. 

Hence volatility/intermittency is a relative concept, and its meaning depends on the 
particular setting under investigation. Once that meaning is clarified the question is how 
to assess the volatility/intermittency empirically and then to describe it in stochastic terms 
and incorporate it in a suitable probabilistic model. 

Such 'additional' random fluctuations generally vary, in time and/or in space, in re- 
gard to Intensity (activity rate and duration) and Amplitude. Typically the volatil- 
ity/intermittency may be further classified into continuous and discrete (i.e. jumps) ele- 
ments, and long and short term effects. 

In finance the investigation of volatility is well developed and many of the procedures 
of probabilistic and statistical analysis applied are similar to those of relevance in turbu- 
lence, for instance in regard to multipower variations, particularly quadratic and bipower 
variations and variation ratios. 

Other important issues concern the modelling of propagating stochastic volatility/inter- 
mittency fields and the question of predictability of volatility/intermittency. 

This paper introduces a concept of realised relative volatility/intermittency and hence 
of realised relative energy dissipation, the ultimate purpose of which is to assess the relative 
volatility/intermittency or energy dissipation in arbitrary subregions of a region C of 
space-time relative to the total volatility/intermittency/energy in C. 

We start the further discussion by describing, in Section 2, how energy dissipation 
in turbulence is defined and traditionally assessed. This is followed by a brief outline 
of some results from the theory of Brownian semistationary (BSS) processes that are 
pertinent for the main topic of the present paper. The definition of realised relative 
volatility /intermittency/energy is given in Section 3. For concreteness and because of 
its particular importance we focus on realised relative energy dissipation. Asymptotic 
probabilistic properties — consistency and a functional central limit theorem — for realised 



For discussions of intermittency and energy dissipation in in turbulence see [Fri95], Chapter 6 and 
[Tsi09], Chapter 7, cf. also the illustration on p. 20 of the latter reference. 
2 Often thought of as Gaussian. 



relative power variations are derived in Section 4. Applications to data on turbulence and 
energy prices are presented in Section 5. Section 6 concludes. 

2 Energy dissipation in turbulence 

In a purely spatial setting the energy dissipation of a homogenous and isotropic turbulent 
field is (up to an ignorable constant involving viscosity) 

1,3=1 •> 

where y\ denotes the velocity at the spatial position x G M 3 . The coarse grained energy 
dissipation over a region C in IR 3 is then given by 

e(C) = ^ j c e(x)dx. 

Furthermore, if only measurements of the velocity component in the main direction x\ of 
the flow are considered one defines the surrogate energy dissipation as 

By Taylor's frozen field hypothesis [Tay38], this may then be reinterpreted as the 
timewise surrogate energy dissipation 

e(t) = {-at 

which would be the relevant quantity in case the measurements were of the same, main, 
component of the velocity but now as a function of time rather than of position. 

Associated to this is the coarse grained energy dissipation corresponding to the interval 
[t, t + u] and given by 

_, v e + {t + u)-e + (t) 
u 
where 

Supposing that the velocity yt has been observed over the interval [0, T] at times 
0, 5, . . . , \T/5\5, when it comes to estimating e + (t), as given by (2), this is traditionally 
done by taking the normalised realised quadratic variation [ys}t/S, where 

\t/8\ 

[ys]t = ^2(yj6-y(j-i)s) 2 - 

3=1 

Correspondingly, the coarse grained energy dissipation over [0, T] is estimated by 

Or) = !^- 

The definitions (1) and (2) of course assume that the sample path y is differentiable. 
On the other hand, going back to Kolmogorov, it is broadly recognised that turbulence can 



only be comprehensibly understood by viewing it as a random phenomenon. 3 Accordingly, 
y should be viewed as a stochastic process, henceforth denoted Y, and it is not realistic 
to assume that its sample paths are differentiable. Thus a broader setting for the analysis 
of the energy dissipation in Y is called for, and in the following we propose and discuss 
such a setting. 

A Brownian semistationary (BSS) process, as introduced in [BNSch09], may be used 
as a model for the timewise development of the velocity at a fixed point in space and in 
the main direction of the flow in a homogeneous and isotropic turbulent field. For focus 
and illustation we shall consider cases where Y is a stationary BSS process, 

Y t = fi+ f g(t- s)a s dB s + f q(t - s)a 2 ds, (3) 

J — oo J— oo 

where g and q are deterministic kernel functions, B is Brownian motion and a is a sta- 
tionary process expressing the volatility/intermittency of the process. In that context the 
gamma form 

git) = ct^e- xt (4) 

of the kernel g has a special role. In particular, if v = 5 and a is square integrable, then 
the autocorrelation function of Y is identical to von Karman's autocorrelation function 
[vKar48] for ideal turbulence. 

In relation to the BSS process (3) with gamma kernel (4) a central question is that 
of determining a 2 from Y. In case the process is a semimartingale the answer is given by 
the quadratic variation of Y; in fact, then [Y] t = o^ + , where 

a 2+ =[a 2 ds (5) 



is the accumulated quadratic volatility over the interval [0, t] . However, in the cases of most 
interest for turbulence, that isf£(^,l)U(l,|) the process Y is not a semimartingale and 
in order to determine a 2+ by a limiting procedure from the realised quadratic variation 

[t/Sj 

\xs\t = J2o r is-Y U . 1)s ) 2 , t>o, 

the latter has to be normalised by a factor depending on 5 and v. Specifically, as shown 
in [BNSch09], this factor is 5c (S)~ where 

c(5) = E{(Gs-G ) 2 } 1/2 (6) 

is defined using the Gaussian core 



G t = [ g(t - s)dB s 



of the process Y. Then have that, as 5 — > 0, 

6 



c{5f 



[Y S ] t A a 2+ . (7) 



Using this result for estimation of a t + requires either that v is known or that a sufficiently 
accurate estimate of v can be found, The latter question has led to detailed studies of the 
application of power and multipower variations to estimation of v, see [BNCP12, BNCP13, 
CHPP13]. 



3 See [Kol41a, Kol41b, Kol41c, Kol62] and, for a recent overview, [Tsi09]. 



3 Realised relative V/I/E 

Supposing again that the velocity Y t has been observed at times 0,5, ... , [T/5\5, we are 
interested in the relative energy dissipation of Y over any subinterval [t, t + u] of [0, T], 
i.e., 

e+(t + u)-e+(t) 

^+(T) ' () 

where £ + (T) is the energy dissipation in [0, T]. Within the turbulence literature, this defi- 
nition of the relative energy dissipation is strongly related to the definition of a multiplier 
in the cascade picture of the transport of energy from large to small scales (see [CSG08] 
and references therein). 

We now introduce the concept of realised relative energy dissipation. Specifically, 
whether Yj is deterministic and differentiable or an arbitrary stochastic process we define 
the realised relative energy dissipation over the subinterval [£, £ + it] as 

Rt(t,t + u) = M]t+u,T-\Y6\t,T> 

where 

[5lt > T " [Y S ] T 

is the realised relative quadratic variation of Y. We note that the quantity R$(t, t + u) is 
entirely empirically based. 

In the "classical" case of turbulence, where Yt is differentiable, as 5 — > we have 

[Y8] t+u -[Ys] t ~5(e + (t + u)-e + (t)) 

and hence, as 5 — > 0, 

RJ(t,t + u)-+ ^^ , 

i.e., the limit equals the relative energy dissipation (8). 

Now suppose that Y is a stationary BSS process (3) with gamma kernel (4) and 
v > \, as it needed for the stochastic integral to exist. Then, if v > \,Y has continuous 
differentiable sample paths, i.e. we are essentially in the "classical" situation. If v = 1 the 
process Y is a semimartingale and the realised quadratic variation [Yg] converges to the 
quadratic variation [Y], that is 

[Y 5 ] t A [Y] t = a 2 t + , 

where a t + is the accumulated quadratic volatility/intermittency (5). Consequently, for 
the realised relative energy dissipation we have 

R+(t, t + u) 4 [Y]t+u " [Y]t - °^ u ~ a " + 



[Y] T 4 + ' 

i.e., the limit is the relative accumulated squared volatility/intermittency. 

Finally, suppose that v G (|, 1) U (1, |), i.e., we are in the non-semimartingale case 
and the sample paths are Holder continuous of order v — 1/2. Then, subject to a mild 
condition on q (see Appendix C for a result covering the case where q is of the gamma 
form), we have again, as 5 — > 0, that 

0-2+ _ rT 2+ 
o+n j- i \ p , t+u °t 

RJ(t,t + u) -)• 2+ 

CTrp 



although 

[Y S ] t « < 



(S-^-^^oo if i/e(|,l) 



^-l) _^ if j/ e (l, f ) 



This follows directly from limiting results in [BNSch09] and [BNCP12]. In view of these 
results, in the turbulence context we view the limit of Rt as the relative energy dissipation. 

Remark 1. As mentioned in Section 2, use of the original assessment procedure (7) requires 
determination of the degree of freedom/smoothness parameter v. The realised relative 
quadratic variation [Y$] t T is entirely empirically determined and its consistency does not 
rely on inference on v. 

4 Probabilistic asymptotic theory of realised relative power 
variations 

We develop now a probabilistic asymptotic theory for realised relative power variations, 
going slightly beyond the earlier discussion of quadratic variations and energy dissipation. 
To highlight the robustness of realised relative power variations to model misspecification, 
we consider both a BSS process and an Ito semimartingale as the underlying process. 
While we limit the discussion to power variations for the sake of simpler exposition, our 
results can be easily extended to multipower variations. 

4.1 Probabilistic setup and consistency 

Let us consider a stochastic process Y = {Y t }t>o, defined on a complete filtered probability 
space (0, T, (J^teM.) P) y i a the decomposition 

Y t = A t + X t , 

where A = {At}t>o is a process that allows for skewness in the distribution of lj. The 
process A is assumed to fulfill one of two negligibility conditions, viz. (10) and (13) given 
below (Appendix C presents more concrete criteria that can be used to check these condi- 
tions). Given a standard Brownian motion B = {B t }tem and a caglad process a = {at}t&.-> 
adapted to the natural filtration of B, we allow for the following two specifications of the 
process X = {X t }t>o- 

(I) X is a local Brownian martingale given by 

X t = I <J s dB s . 
Jo 

(II) X is a BSS process given by 

X t = g(t - s)a s dB s , 

J — oo 

where g : (0, oo) — >• R is a square integrable weight function such that 

/g(t — s) o~ s ds < oo 
for all t > 0. 



a.s. 



On the one hand, by choosing A to be absolutely continuous in the case (I), we see that 
this framework includes rather general Brownian semimartingales. On the other hand, in 
the case (II) the process Y is typically a non-semimartingale, as discussed above. 

Recall that for any p > 0, the p-th. order realised power variation of the process Y with 

lag 5 > is given by 

[t/s\ 

3=1 

The power variations [As] an d [-Xa]^' are, of course, defined analogously. Similarly 
to earlier literature on power and multipower variations of BSS processes, [BNCP12, 
BNCP13, CHPP13], we assume that the kernel function g behaves like f near zero for 
some i/£(fj,1)U(1,|), or more precisely that 

g{t)=t»- 1 L g (t), 

where L g is slowly varying at zero, which implies that X is not a semimartingale in the 
case (II). Then, under some further regularity conditions 1 on g, which are satisfied when 
g is the gamma kernel (4), and defining the normalisation factor c(S) by (6), we have 

-^rai^^r, (9) 

where a P+ = jj \a s \ p ds and m p = E{|£| p } for £ ~ JV(0, 1), by Theorem 3.1 of [CHPP13]. 
In the case (I), setting e(<5) = y/5, the convergence (9) holds without any additional 
assumptions, e.g., by Theorem 2.2 of [BGJPS06]. Additionally, note that the convergence 
(9) holds also when X is replaced with Y, provided that 

S -[A s ]?-^0. (10) 



c{5)p 1 on s^o 

For fixed time horizon T > 0, we introduce the p-th order realised relative power 
variation process over [0, T] by 

\y7] {p) - [ ^- 0<t<T 

The relative power variation has the following evident consistency property. 
Theorem 2. If (9) and (10) hold, then 

uniformly in t G [0, T]. 

Remark 3. The relative integrated volatility/intermittency a p T can be seen as a "cumula- 
tive distribution function" of volatility/intermittency on [0, T]; we have Oq T = 0, Oj, T = 1, 
and 1 1—)- a P T is non-decreasing. 

Remark 4. Theorem 2 can be generalized to a spatial setting, applying specifically to 
ambit fields driven by white noise, using the asymptotic results in [Pakl3]. 



4 See [CHPP13]. Most importantly, they include the assumption that c(8) = <5" 1//2 L C (<5), where L c is 
slowly varying at zero. 



4.2 Central limit theorem 

Recall first that random elements XJ\ , U% , . . . in some metric space U converge stably (in 
law) to a random element U in U, denned on an extension (fi', J-',P') of the underlying 
probability space (Q, J 7 , P), if 

W(Un)V} ► E'{f(U)V} 

n— >-oo 

for any bounded, continuous function / : U — > M. and bounded random variable V on 
(O, J-, P). We denote stable convergence by — >. Clearly, stable convergence implies or- 
dinary convergence in law, while the converse is not true. This mode of convergence was 
introduced by Renyi [Ren63]. 

Remark 5. The usefulness of stable convergence can be illustrated by the following example 

st 

that is pertinent to the asymptotic results below. Suppose that U n — > OS, in M. as n — > oo, 
where £ ~ N(0, 1) and 9 is a positive random variable independent of £. In other words, 
U n follows asymptotically a mixed Gaussian law with mean zero and conditional variance 
8 2 . If 9 n is a positive, consistent estimator of 9, i.e., 9 n — > as n — > oo, then the stable 

convergence of U n allows us to deduce that U n /9 n — > iV(0, 1). We refer to [Ren63, AE78] 
and [JS03, pp. 512-518] for more information on stable convergence. 

Let us write D([0, T]) for the space of cadlag functions from [0, T] to R, endowed with 
the usual Skorohod metric [JS03, Chapter V]. (Recall, however, that convergence to a 
continuous function in this metric is equivalent to uniform convergence.) Under slightly 
strengthened assumptions, the realised power variation of X satisfies a stable central limit 
theorem of the form 

rl/2 (# W ' (P) "' vff ) 7f^ V^j[V.|W, mD([0,T]) ; (11) 

where \x,p > is a deterministic constant and {Wjjtgro^] a standard Brownian motion, 
independent of J 7 , defined on an extension of (Q, J 7 , P). Indeed, in the case (I) the con- 
vergence (11) holds with \x,p = w-2p — mz, provided that a is an ltd semimartingale (see 
Theorem 2.4 of [BGJPS06]). Moreover, we have (11) also in the case (II) if we make 
the restriction v e (o> 1) — the situation of most interest concerning turbulence — and 
assume that a satisfies a Holder condition in expectation (see Theorem 3.2 of [CHPP13]). 
Then, in contrast to the semimartingale case, 

Ax, P = Ap(z^) > m 2p - ml, (12) 

where X p : (|, 1) — > (0, oo) is a continuous function defined using the correlation structure 
of fractional Brownian noise (see Appendix B for the definition and proof of continuity). 
Analogously to (10), the convergence (11) extends to the power variation of Y if 

V^ r , l(p) 



c W ^l' t?°- (13) 

The realised relative power variation of Y satisfies the following central limit theorem, 
which is an immediate consequence of Lemma 10 in Appendix A. 



Theorem 6. If (11) and (13) hold, then 



8^0 



m 



Jo l^l Pds 



V Jo 



[Ys]tT-^T )t^ 7T7^ [ \<T a \ p dW a -dft \cs\ p dW s \ (14) 



T 



in D([0,T]), where Xx, P and W are as in (11). 



Remark 7. In the case (II) the restriction v G (|, 1) can be relaxed when one considers 
power variations defined using second or higher order differences of Y (cf. [BNCP13, 
CHPP13]). Then, (11) holds for all v G (|,1) U (1, §). As Theorems 2 and 6 do not 
depend on the type of differences used in the power variation, they obviously apply also 
in this case. 

Conditional on J 7 , the limiting process on the right-hand side of (14) is a Gaussian 
bridge. In particular, its (unconditional) marginal law at time t G [0, T] is mixed Gaussian 
with mean zero and conditional variance 



{m v a, 



%- ( (1 - o^fo-r + Kt) 2 (4 P+ - oF)) • (15) 



Note also that when a is constant, the limiting process reduces to a Brownian bridge. In 
effect, the result is analogous to Donsker's theorem for empirical cumulative distribution 
functions (see, e.g., [Kos08] for an overview of such results). 

Clearly, we may estimate the asymptotic variance (15) consistently using 

v,(s) = - — ^-^((i - mtl) 2 nt' ] + (Si2) 2 ([<» - wiS 2 *')). 

8 -mo„ • \\Ya\x,') \ x 7 K y J 



8-m 2p - ([Y S ] T 

In the case (II) the estimator Vt(8) is not feasible as such since v appears as a nuisance 
parameter in \x,p = \p(v)- However, we may replace Xx, P with X p (i>s), where vs is any 
consistent estimator of v based on the observations Yq, Y$, . . . ,Y\t/s\6 (t ne Y have been 
developed in [BNCP12, BNCP13, CHPP13]). Using the properties of stable convergence, 
we obtain the following feasible central limit theorem. 



Proposition 8. If (11) and (13) hold, then for any t G (0, T), 

->iV(0,l). 



s- 1/2 ([ys]tl-^Ti „ 



VW) 



8^0 



4.3 Inference on realised relative V/I/E 

Proposition 8 can be used to construct approximative, pointwise confidence intervals for 
the relative volatility /intermittency &t T . Since, by construction, a^ T assumes values in 
[0, 1], it is reasonable to constrain the confidence interval to be a subset of [0, 1]. Thus, 
we define for any a G (0, 1) the corresponding (1 — a) ■ 100% confidence interval as 

m3x^[Y s ] t>T - Z!_ a / 2 ■ y/8Vt(8),0J, mm{[Y 5 ] t;r + Zl _ a/2 ■ y/8V t (8),lj , 

where Zi- a /2 > is the 1 — |-quantile of the standard Gaussian distribution. 

Another application of the central limit theory is a non-parametric homoskedasticity 
test that is similar in nature to Kolmogorov-Smirnov and Cramer-von Mises goodness- 
of-fit tests for empirical distribution functions. This extends the homoskedasticity tests 
proposed by Dette, Podolskij and Vetter [DPV06] and Dette and Podolskij [DP08] to a 
non-semimartingale setting. To formulate our test, we introduce the hypotheses 

J H :a t = o- for all t G [0,T], 
1 Hi : at 7^ 0"o for some t G [0, T\. 



Theorem 6 implies that under Hq, 



It is intuitively appealing to use the L 1 -norm to assess the distance of realised relative 
power variations — although other metrics or norms could also be used. Thus, we define 
the test statistic as 

dp) t 



rT 

m, 



[Ys]t,T - f 



dt. (17) 



y/6 ■ T ■ X x ,p 

Also in (17), we may use A p (i> < j) instead of Xx :P in the case (II). By (16) and the scaling 
properties of Brownian motion, we have under Ho, 

where {W^jjgjo.i] is a standard Brownian bridge independent of J 7 , while under Hi the 
statistic S$ diverges to infinity. The cumulative distribution function of the functional on 
the right-hand side of (18) has been derived by Johnson and Killeen [JK83] and it can be 
stated using the Airy function, which allows for a straightforward numerical implementa- 
tion. 

5 Applications to empirical data 

5.1 Brookhaven turbulence data 

We apply the methodology developed above first to data of turbulence. The data consist 
of a time series of the main component of a turbulent velocity vector, measured at a 
fixed position in the atmospheric boundary layer using a hotwire anemometer, during an 
approximately 66 minutes long observation period at sampling frequency of 5 kHz (i.e. 
5000 observations per second). The measurements were made at Brookhaven National 
Laboratory (Long Island, NY), and a comprehensive account of the data can be found in 
[DrhOO]. 

As a first illustration, we study the observations up to time horizon T = 800 millisec- 
onds. Using the smallest possible lag, 5 = 0.2 ms, this amounts to 4000 observations. 
Figure 1(a) displays the squared increments corresponding to these observations. As 
a comparison, the same time horizon is captured in Figure 1(b) but with lag 5 = 0.8 
ms. Figure 1(c) compares the associated accumulated realised relative energy dissipa- 
tions/quadratic variations. The graphs for these two lags show very similar behaviour, 
exhibiting how the total time interval is divided into a sequence of intervals over which 
the slope of the energy dissipation is roughly constant. On the other hand, the amplitudes 
of the volatility/intermittency are of the same order in the whole observation interval. 

To be able to draw inference on relative volatility/intermittency using the data, we 
need to address two issues. Firstly, for this time series, the lags 5 = 0.2 ms and 5 = 0.8 ms 
are below the so-called inertial range of turbulence, where a BSS process with a gamma 
kernel --a model of ideal turbulence -- provides an accurate description of the data 
(see [CHPP13], where the same data are analysed). Secondly, the data were digitised 
using a 12-bit analog-to-digital converter. Thus, the measurements can assume at most 
2 12 = 4096 different values, and due to the resulting discretisation error, a non-negligible 
amount of the increments are in fact equal to zero (roughly 20 % of all increments) . These 

10 
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Figure 1: Left panel, Brookhaven turbulence data: (a) The squared increment process 
with lag S = 0.2 ms over the time horizon T = 800 ms. (b) The squared increment process 
with lag 5 = 0.8 ms over the same time horizon T = 800 ms. (c) The realised relative 
energy dissipation/quadratic variations corresponding to 5 = 0.2 ms and 5 = 0.8 ms, and 
the same time horizon, T = 800 ms, as in plots (a) and (b). Right panel, logarithmic 
EEX electricity spot prices: (d) The squared increment process with lag 5 = 1 day over 
the time horizon T = 1775 days, (e) The squared increment process with lag 5 = 4 days 
over the same time horizon T = 1775 days, (f) The realised relative quadratic variation 
corresponding to 5 = 1 day and 5 = 4 days and the same time horizon, T = 1775 days, as 
in plots (d) and (e). 
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Figure 2: Brookhaven turbulence data: (a) Estimates of u, using the change-of-frequency 
method and lag 5 = 80 ms, for each one-minute subperiod and the value predicted by 
Kolmogorov's (K41) scaling law. (b) Test statistics and critical values for the constancy 
of a for each subperiod. The red bars indicate the 40th and 44th subperiods, which are 
analysed in more detail in Figure 3. 

discretisation errors are bound to bias the estimation of the parameter v, which is needed 
for the inference methods. We mitigate these issues by subsampling, namely, we apply 
the inference methods using a considerably longer lag, 5 = 80 ms, which is near the lower 
bound of the inertial range for this time series [CHPP13, Figure 1]. 

We divide the time series into 66 non-overlapping one- minute- long subperiods, testing 
the constancy of a, i.e., the null hypothesis Ho, within each subperiod. Figure 2(a) displays 
the estimates of v for each subperiod using the change-of-frequency method [BNCP13, 
CHPP13]. All of the estimates belong to the interval (^ , 1) and they are scattered around 
the value v = I predicted by Kolmogorov's (K41) scaling law of turbulence [Kol41a, 
Kol41c]. The homoskedasticity test statistics and critical values, derived using (18), in 
Figure 2(b) indicate that the null hypothesis of constancy of a is typically rejected. To 
understand what kind of intermittency the test is detecting, we look into the extremal 
cases, the 40th and 44th subperiods (red bars). To this end, we plot the realised relative 
energy dissipations, with 5 = 80 ms, during the 40th and 44th subperiods in Figure 3(a) 
and (b), respectively. We also include the pointwise confidence intervals, the p- values 
of the homoskedasticity test, and as a reference, the realised relative energy dissipations 
using the smallest possible lag 5 = 0.2 ms. While the realised relative energy dissipations 
exhibit a slight discrepancy between the lags 5 = 80 ms and 5 = 0.2 ms, it is clear that 
40th subperiod contains significant intermittency, whereas the during the 44th subperiod, 
the (accumulated) realised relative energy dissipation grows nearly linearly, apart from 
the final 10 seconds. 
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Figure 3: Brookhaven turbulence data: Realised relative energy dissipation during the 
40th (a) and 44th (b) subperiods with 5 = 80 ms and 5 = 0.2 ms. Additionally, p-values 
for the hypothesis Ho, estimates of v using the change-of- frequency method, and 95% 
pointwise confidence intervals, all using the lag 5 = 80 ms. 

5.2 EEX electricity spot prices 

We also briefly exemplify the concept of relative volatility using electricity spot price data 
from the European Energy Exchange (EEX) . Specifically, we consider deseasonalised daily 
Phelix peak load data (that is, the daily averages of the hourly spot prices of electricity 
delivered between 8 AM and 8 pm) with delivery days ranging from January 1, 2002 to 
October 21, 2008. Weekends are not included in the peak load data, and in total we have 
1775 observations. This time series was studied in [BNBV12] and the deseasonalisation 
method is explained therein. As usual, we consider here logarithmic prices. 

Figure 1(d) shows the squared increments up to the total time horizon T = 1775 days 
with lag 5 = 1 day. The same time horizon is captured in Figure 1(e) but with a reso- 
lution 5 = 4 days. Figure 1(f) compares the corresponding accumulated realised relative 
quadratic variations. The results for these two lags do not show the same similarity as 
with the turbulence data (Figure l(a-b)). Judging by eye, we observe that the intensity 
of the volatility is changing with lag 5. This lag dependence is also observed in the ampli- 
tudes, again in contrast to the figures on the left hand side. (However, more quantitative 
investigation of such amplitude/density arguments is outside the scope of the present pa- 
per.) The dependence of the estimation results on the lag 5 is, at least partly, explained 
by the relatively low sampling frequency of the data. With 5 = 1 day, the increments are 
dominated by a few exceptional observations (which may correspond to jumps or intraday 
volatility bursts). Choosing 5 = 4 days reduces the contribution of these observations 
since the time series exhibits significant first-order autocorrelation [BNBV12, Figure 1]. 

Remark 9. It was shown in [BNBV12] that by suitably choosing both g and q to be of 
gamma type it is possible to construct a BSS process with normal inverse Gaussian one- 
dimensional marginal law, which corresponds closely to the empirics for the time series of 
log spot prices considered. Moreover, the estimated value of the smoothness parameter v 
for this time series falls in the interval (3, 1). 
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6 Conclusion 

The definition of realised relative energy dissipation introduced in this paper applies to 
any continuous time, real valued process Y. An extension to vector valued processes is an 
issue of interest, in particular in relation to the definition (1) of the energy dissipation in 
three-dimensional turbulent fields. 

The extent to which the realised volatility/intermittency/energy is an empirical coun- 
terpart of what can be conceived theoretically as relative volatility/intermittency/energy 
depends on the model under consideration. As discussed above this is the case, in particu- 
lar, both under Brownian semimartingales, as these occur widely in mathematical finance 
and financial econometrics, and under stationary BSS processes. 

In the timewise stationary setting, the realised relative energy dissipation is a param- 
eter free statistic which provides estimates of the relative energy in subintervals of the 
full observation range, by relating the quadratic variation over each subinterval to the 
total realised energy for the entire range. It provides robust estimates of the relative ac- 
cumulated energy as this develops over time and is intimately connected to the concepts 
of volatility/intermittency and energy dissipation as these occur in statistical turbulence 
and in finance. This was illustrated in connection to the class of BSS processes with g of 
the gamma form. 
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A Stable convergence of normalised processes 

Theorem 6 is a corollary of the following simple result concerning the stable convergence 
of a process that has been normalised by its terminal value. 

Lemma 10. Let X,X l ,X 2 ,... be random elements in D([0,T]), with non- decreasing 
sample paths, defined on (Q, J 7 , P) and let £ be a random element in D([0,T]) defined on 
an extension (f2', J-',P') of (£l,J-',P). Suppose, moreover, that almost surely S, G C([0, T]), 
X G C([0, T]), X£ ^ for any n G N, and X T ^ 0. If 

n(x? - X t ) -±+ & inD([0,T]), (19) 



then 



-ni^-^)r^^((,-^ir) » D ([0,T,). 



x% x T/ 

Proof. Since X n and X have non-decreasing paths and X G C([0,T]), we have 

X?_Xt_ 

\Xt\ o / /' 



sup 

0<t<T 



xz X T 



< ^^ sup |A7 - X t \ — £->• 



by (19). Due to the properties of stable convergence, we obtain 

X t\ st ( t X t \ , r , 2 



S(AT-X,),^j^^,^j m O ([0,rr. (20) 
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Let us now consider the decomposition 

Using again the fact that convergence to a continuous function in D([0,T]) is equivalent 
to uniform convergence, it follows that the map (x,y) i— y x — x(T)y from D([0,T]) 2 to 
D([0,T]) is continuous on C([0,T]) 2 . Since f G C([0,T]) and X G C([0,T]) a.s., the 
assertion follows from (20) and the properties of stable convergence. □ 

B Continuity of the correlation structure of fractional Brow- 
nian noise 

The function X p : (|, 1) — > (0, oo) that appears in (12) is defined by 

OO i OO \ 

1=2 V j=l ' 

where 02, 03, . . . are the coefficients in the expansion of the function u p (x) = \x\ p — m P , 
x £ R, in second and higher-order Hermite polynomials x 2 — 1, x 3 — 3x, . . . , satisfying 
Yl'i=2 ^■ a l < °°- I n the mos t important case p = 2 we have, clearly, 02 = 1 and a; = 
for all Z > 2. The sequence (p v (j))?^=i coincides with the correlation function of fractional 
Brownian noise with Hurst parameter v — 1/2, namely 

Pu(j) = \{{j + I) 2 "" 1 - 2i 2 "~ 1 + (j - l) 2 ^ 1 ), j > 1. 

The following lemma justifies the use of \p{v§) as a consistent estimator of \x,p in the 
case (II). 

Lemma 11. The function \ p is continuous. 

Proof. Let us first show that for any I > 2, the map 

00 

"-►EM;)' (21) 

is continuous on (^l)- Note that v 1— >• /5^(j) is clearly continuous for any j > 1 and, 
moreover, that |pj/(j)l' ^ Pv{j) 2 since |/v(j)| ^ 1. Following [DM03, p. 419], write 

Mi) = \j 2u - l 9u(l/j), (22) 

where ^(x) = (1 + x) 2u ~ l — 2 + (1 — x) 2u ~ l , x G R. A straightforward computation shows 
that s^(0) = 0, #^(0) = 0, and 

<^'(x) = (21/ - 2) (2!/ - 1)((1 + x) 2 ^ 3 + (1 - x) 2 ^ 3 ). 

In particular, we may define C = sup„ e n / 2) i) sup a;e fo ) i/2] |<?"(x)| < 00, whence |^(rc)| < 
Cx 2 for all v G (i, 1) and x G [0, ^]. Applying this bound to (22) implies that 

Mj)l < Cj 2 - 3 < -, j > 2, 
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for any v € (i, 1). The continuity of the function in (21) follows now from Lebesgue's 
dominated convergence theorem. Finally, the actual assertion follows by invoking the 
bounds 



OO / OO 1 v 

<Ei^)I 2 <^(i + e4)< 

3=1 V 3=2 J J 



OO, 



3=1 

which hold for alH > 2, and Lebesgue's dominated convergence theorem again. □ 



C Sufficient conditions for the negligibility of the skewness 
term 

Suppose first that the process A = {A t }t>o is given by 

A t = [i + / a s ds, 
Jo 

where \i G R is a constant and the process {at}t>o is measurable and locally bounded. 
Then we can establish rather simple conditions for its negligibility in the asymptotic results 
for power variations. By Jensen's inequality, we have for any p > 1, s > 0, and t > 0, 

\A s -A t \P<C a -\s-t\P, 

where C a > is a random variable that depends locally on the path of a. Thus, the 
condition (10) holds whenever 

W)^ ' <23) 

and (13) holds if 

AP-l/2 

Notably, in the semimartingale case (I), i.e. c(5) = v6, (23) is always true and (24) holds 
if p > 1. 

Suppose now, instead, that A follows 

A t = n+ q(t- s)a s ds, (25) 

J — OO 

where q is the gamma kernel 

q (t) = c't^- l e- pt 

for some c' > 0, rj > 0, and p > 0. We assume that the process {ai}i g K is measurable, 
locally bounded, and satisfies 

/u 
q(u — s)\a s \ds < oo a.s. (26) 

-OO 

for any t > 0, which is true, e.g., when the auxiliary process J_ q(u — s)\a s \ds, u > 0, 
has a cadlag or continuous modification. 

Lemma 12. Suppose that A is given by (25) and that (26) holds. Then in the case (II) 
i/ie condition (10) ZioWs i/min{?7, 1} > v — \ and (13) holds i/min{n, 1} > z/ — ^— . T/ie 
corresponding restrictions in the case (I) can be obtained by setting v = 1 above. 
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Proof. Let us first look into the properties of q. For the sake of simpler notation, we make 
the innocuous assumption that d = 1. Since 

«'(*)= (^-p)?(*), (27) 

we find that q is decreasing when r\ < 1. When 77 > 1, g is increasing on (0, 2= — ) and 

decreasing on (2^-, 00). 

Let £ > be fixed, 5 G (0, 1), and let j > 1 be such that jS < t. Below, all big O 
estimates hold uniformly in such j. We consider the decomposition 

HS r(j-l)S 

Ajs-Af j _ 1 )$= q(j5 - s)a s ds + (q(j5 - s) - q((j - 1)5 - s))a s ds 

J {3-1)6 J(j-2)6 

rU-ty 
+ / {q(jS ~ s) - q({j - 1)5 - s))a s ds 



+ / {q(j6 -a)- q((j - 1)5 - s))a s ds 
Il+Ij+Il + Il 



max < , 1 



where 



P 

When 77 > 1, q is bounded and we have \ll + J? | = a^O(5), where 

a* t = sup \a s \ < 00 a.s., 

s*<s<t 

and when r/ < 1, we find that 

r 6 



I}+I]\<2a* t f q(s)ds = a* t O(5V). 
Jo 



/o 
Next, we want to show that 

|J|| = ajO(<5 min ^' 1 >). (28) 

In the case 77 > 2 the derivative g' is bounded and (28) is immediate. Suppose that ij < 2. 
Then, |g'(i)| < Ct v ~ 2 on any finite interval, where C > depends on the interval. Using 
the mean value theorem, we obtain 

\ll\<Ca$6 {(j-l)5-s) v - 2 ds, 

J s* 

which implies (28). To bound |l||, note that, by (27), \q'(t)\ < C'q(t) for all t > -s*, 



■n-i 

jiiOLcxiilj. r ui a-J-iJ 0^0, vvc ncxvc ^J — x ju — o ^ - 

mean value theorem 



where C > is a constant. For any s < s* , we have (j — 1)5 — s > - — . Thus, by the 



I {q(j6 -s)- q((j - 1)5 -s))\< C'q((j - 1)5 - s)5 

and, consequently, 

-0-1)5 



\ll\<C'5J' J q{(j-l)5-s)\a s \ds = A* t O(5). 
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Collecting the estimates, we have 

\A jS - A u _ 1)5 \ = maxK,^}0(5 min{ "' 1} ) 
uniformly in j, whence 

Checking the sufficiency of the asserted conditions is now a straightforward task. 



□ 
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